arXiv:1506.04960v2 [hep-ex] 13 Nov 2015 


Chinese Physics C Vol. XX, No. X (201X) XXXXXX 


Search for '?/i(45') —> r/J/?/; in ^ rjJ/tljK^ and e“^e processes 

Gao Xu-yang^ Wang Xiao-long^’® Shen Cheng-ping^’^' 

^ School of Physics and Nuclear Energy Engineering, Beihang University 
Xueyuan Road No.37, Haidian District, Beijing 100191, China 
^ Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 
®CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 

Abstract: We search for the '!^(4S') state in the l rjJ/ipK^ and e"'"e“ vJ/'4’ processes based on the Belle 

measurements with the assumed mass M = (4230 ±8) MeV/c^ and width P = (38 ±12) MeV. No significant signal 
is observed in the tjJ/V’ mass spectra. The 90% confidence level upper limit on the product branching fraction 
B{B^ — > riJ/ip) < 6.8 x 10“® is obtained in the B^ —l rjJ/ipK^ decays. By assuming the partial 

width of ip{4:S) —>■ e^e“ to be 0.63 keV, a branching fraction limit B{ip{4S) —>■ rjj/ip) < 1.3% is obtained at the 90% 
confidence level in e'''e~ -^r]J/ip, which is consistent with the theoretical prediction. 
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1 Introduction 

Potential models predict some charmonium states 
above the DD threshold, but the number of observed 
states in experiments is more than the number predicted. 
The states which have been seen outwith the theoretical 
predictions are normally referred to as exotic states or 
XYZ particles. Many XYZ states have been announced 
in various processes, for example, the observation of 
Y(3872) in B decays fl, the Y(4260) 0, Y(4360) 0 and 
Y(4660) 0 in e+e“ annihilation, and the Y(3915) [sj ob¬ 
served in the two-photon process. 

On the other hand, there are still some charmonium 
states predicted by the potential models which have not 
yet been observed experimentally, especially in the mass 
region higher than 4 GeV/c^, such as r]c{4S), 

ip{4:S) and ip{5S). To some degree, some XYZ states 
are regarded as candidates for these unfound predicted 
states. 

Searching for these missing predicted states is very 
helpful to test the potential models. When checking the 
mass spectra of the observed charmonia with spin-parity 
= 1 and comparing them with those of the cor¬ 
responding bottomonia, there might be a charmonium 
state 54 ( 45 ') at about 4.2 GeY ^ compared to the T(45) 
state [6| . The authors in Ref. [Gj predicted that this miss¬ 
ing charmonium state has a mass of 4.263 GeV /(? and a 
very narrow width. As a state with the same spin-parity 
1 , the Y(4220) 0 may be a good candidate for the 


'tp{4:S) state. 

Recently, the BESIII Gollaboration performed a 
study on the decay e+e“ —>■ (J = 0,1,2) 0, where 

the Born cross sections at 9 energy points were mea¬ 
sured. When using a Breit-Wigner (BW) function to 
fit the experimental data of e+e“ —>■ ujXco, a resonant 
structure with mass M = (4230 ±8) MeV/c^ and width 
r= (38±12) MeV was observed with a statistical signif¬ 
icance more than 9cr. However, for the remaining pro¬ 
cesses e+e“ — >■ oJXci and e+e“ —>■ wXc2j there were no 
significant signals. 

To understand this novel phenomenon, different ex¬ 
planations of this resonance were given, which included 
a tetraquark state 0, the missing higher charmonium 
state V’(45) 0, and the known charmonium resonance 
4(4160) [i3|. 

The authors in Ref. 0 checked the thresholds of 
uJXco, ^Xci and uJXc 2 , which are 4.197 GeV/c^, 4.293 
GeV/c^ and 4.338 GeV/c^, respectively. The central 
mass of 4(45) is just above the wXco threshold and below 
the u;Xci ,2 thresholds. Accordingly the newly observed 
structure in e"* *'e“ —>■ wXco could be the missing charmo¬ 
nium 4(45) state, and the e+e“ —>■ a;Xci ,2 processes are 
kinematically forbidden to 4(45). Stimulated by this, 
the authors estimated the meson loop contribution to 
4(45) —>■ uJXco and found the evaluation can overlap with 
the experimental data in a reasonable parameter range. 
As a typical transition accessible by experiment, the de¬ 
cay 4(45) —>■ rjJ/ip similar to ip{^S) —>■ wXco can occur. 
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So the authors in Ref. @ also extended the theoreti¬ 
cal calcution to ^/’(4S') —?> ryJ/'f/' and predicted the upper 
limit on the branching fraction of ^/’(4S') —?> TyJ/'i/i to be 
less than 1.9 x 10“® via the hadronic loop mechanism [^. 

As indicated in Ref. Q, the predicted upper limit 
of '0(45') —>■ 77J/0 can be accessible at Belle and the 
forthcoming Bellell. We noticed that the Belle experi¬ 
ment previously measured the > rjJ/tjjK^ [H and 
e'''e“ —>■ jyJ/'i/’ [l2| processes, where the /yJ/0 invariant 
mass distributions were given. Hence, in this work we 
fit the rjJ/ip mass spectra from the —?> rjJ/ipK^ and 

e+e“ —>■ r]J/'ip processes to search for the 0(45') state. 
The experimental measurements can be taken as a test 
of the theoretical calculation. 

This work is organized as follows. We present the 
detailed fit results to the ?y J/0 mass spectra from B^ —^ 
7yJ/0Ar=*= and e+e“ —>■ ryJ/0 processes with the 0(45) 
state included in Sec. 2 and Sec. 3. If no clear 0(45) 
signal is observed, the branching fraction limits at the 
90% confidence level (C.L.) will be given with the sys¬ 
tematic errors included. The last section ends with the 
conclusion and discussion. 


value of the likelihood as a function of Wig- To take 
into account the systematic uncertainty, the above likeli¬ 
hood is convolved with a Gaussian function whose width 
equals the total systematic uncertainty described below. 
The upper limit on the number of 0(45) signal events is 
22.7 at 90% C.L. 



Fig. 1. (color online) The r]J/'ip invariant mass dis¬ 
tribution from -qj/ipK^ decays. The dots 

with error bars are from data, the solid curve is 
the best fit for the total signal and the dotted line 
shows the fitted background shape. 


2 Search for 0 ( 45 ') in B Decays 


Using 772 X 10® BB pairs collected with the Belle de¬ 
tector, the decays B^ —^ qJ/ipK^ were studied to search 
for a new narrow charmonium(-like) state X in the 77 J /0 
mass spectrum, where the J/0 and q mesons were re¬ 
constructed by a lepton-pair (I = e, [x) and two 

photons 11|. Except for the known 0' —>■ ryJ/0 decay, 
no significant narrow excess was found in the qJ/'tp rnass 
spectrum. 

Figure [T] shows the qJ/ip mass distribution of interest 
after all the event selection requirements are applied. A 
binned maximum likelihood fit to the ryJ/0 mass distri¬ 
bution is performed to extract the signal and background 
yields. A BW function (mass and width fixed at 4.23 
GeV/c^ and 38 MeV i) is convolved with a Gaussian 
function (the mass resolution is about 11 MeV/c^) as the 
0(45) signal shape and a second polynomial function is 
taken as the background shape. The fit range and results 
to the qJ/ij) mass spectrum are shown in Fig. [1] 

From the fit, we obtain 5.9±5.5 signal events, with a 
statistical significance of 0.9a, from the difference of the 
logarithmic likelihoods, —21n(£o/^max), taking the dif¬ 
ference in the number of degrees of freedom (Andf = 1) 
in the fits into account, where £0 and C^ax are the like¬ 
lihoods of the fits without and with a resonance compo¬ 
nent, respectively. 

We determine a Bayesian 90% C.L. upper limit on the 
number of 0(45) signal events (A^sig) by finding the value 

fVjT such that CdN^ig = 0.90, where 

fVsig is the number of 0(45) signal events and C is the 


There are several sources of systematic error for the 
branching fraction measurement. Most of the system¬ 
atic errors are the same as those in Ref. [H , except that 
the dominant uncertainty associated with the fitting pro¬ 
cedure is different, which is estimated by changing the 
order of the background polynomial, the range of the 
fit, the 0(45) mass and width by ±lcr. Finally, the un¬ 
certainty due to the fitting procedure is 11%. Assum¬ 
ing all the sources are independent and adding them in 
quadrature, the final total systematic uncertainties are 
summarized in Table [TJ 

Table 1. Relative systematic errors (%) on the 
product of the branching fraction B{B^ —>■ 
0(45)A±)5(0(45) 7 ?J/0). 


Source 

Relative error (%) 

Tracking efficiency 

1.1 

Lepton identification 

2.4 

Charged kaon identification 

1.4 

T) —>• 77 efficiency 

3.0 

Signal MC simulation statistics 

0.5 

Secondary B 

0.7 

^BB 

1.4 

Fitting procedure 

11 

Total 

12 


The 90% C.L. upper limit is set on the product 
branching fraction B{B^ —0(45)Ar=*=)i3(0(45) —>■ qJ/ix) 
using 


Nbb X e X ^ i+£~) X B{q^^"f) ’ 
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where , Ngg, e, ^ £'^£~) and are 

the upper limit on the number of ^/’(4S') signal events at 
90% C.L., the number of BB pairs, the corrected detec¬ 
tion efficiency of 9.23% at 4.23 GeV/c^ obtained from the 
fitted efficiency curve using the efficiencies at ' 0 ', i/'CdOdO) 
and '0(4160) points 11[, the branc hing fractions of J/0 
to lepton pair and 77 to two photons [l^, respectively. Fi¬ 
nally, the 90% C.L. upper limit on the product branching 
fraction B{B^ — 7 > 0(4S')iC=*=)S(0(4S') —>■ r]J/'ijj) is found 
to be 6.8 X 10“®. 


3 Search for '^(45') in e+e -^rjj/ip 


The cross section for e'*'e“ —^ rjj/ip between ^/s = 3.8 
and 5.3 GeV was measured using 980 fb“^ of Belle data, 
where the 77 was reconstructed with its 77 and 
decays and J/0 was reconstructed via its leptonic de¬ 
cays. Two distinct resonant structures, the 0(4040) and 
0(4160), were observed [l^ . 

To obtain the transition rates of 0(4040) and 0(4160) 
to the 77 J/0 final state, an unbinned maximum likelihood 
fit was performed to the 77 J/0 mass spectra from the sig¬ 
nal candidate events and the 77 and J/0 sideband events 
simultaneously 12| . The fit to the signal events includes 
two coherent P-wave BW functions convolved by the ef¬ 
fective luminosity and efficiency curve for 0(4040) and 
0(4160) signals and an incoherent second-order poly¬ 
nomial background; the fit to the sideband events in¬ 
cludes the same background function only. Due to the 
low statistics, the masses and widths of the 0(4040) and 
0(4160) were fixed 1^ and the effects of mass resolution 
were small and therefore were neglected [l^ . 

Similarly here, to obtain the transition rate of 0(45*) 
to 77 J /0 final state, a binned maximum likelihood fit 
with three coherent P-wave BW functions for 0(4040), 
0(4160) and 0(45') is applied to the e+e“ cross 

sections directly, as shown in Fig. [2j In the fits, be¬ 
sides the masses and widths of 0(4040) and 0(4160) 
are fixed [ill, the 0(45) parameters are also fixed 
Figure [2] and Table [2] show the fit results. There are 
four solutions with equally good fit quality. The results 
of PFe+e- for 0(4040) and 0(4160) are consistent with 
the published results within errors [l2j| . The significance 
of the 0(45) is estimated by comparing the likelihood 
of fits with and without 0(45) included. We obtain 
a statistical significance of 2.6(T. The most conserva¬ 
tive upper limit with the systematic errors included on 
P(0(45) —t 77 J/ 0 )F^|'^®^ is obtained to be 8.2 eV at 90% 
C.L., which corresponds to Solution III in Fig.[2Dc). 

Most of the systematic errors in the 13(0(45) — 
77 J/ 0 )F^ i'^^'^^ measurement are the same as those in 
Ref. [l2| except the dominant systematic error from fit 
uncertainty (64%), which incfudes the uncertainties on 
the mass and width of 0(45) state by changing the nom¬ 


inal values by Icr Q, and the fit range. Assuming all the 
sources are independent and adding them in quadrature, 
we obtain total systematic error in PFe+e- of 65% for 
0(45), as shown in Table [31 



Fig. 2. (color online) The e'’'e“ -^rjj/xj} cross sec¬ 
tion distributions and the fit results described in 
the text. There are four solutions with the co¬ 
herent 0(4040), 0(4160) and 0(45) signals. The 
curves show the best fit for the measured cross 
sections and the contribution from each BW com¬ 
ponent. The interference term for each solution is 
not shown. 

Table 3. Relative systematic errors (in %) in the 
5(0(45)—)■ 77 J/0)r)f/:‘^/.^ measurement. 


Source 

Relative error (%) 

Particle identification 

5.5 

Tracking efficiency 

1.4 

Photon reconstruction 

4.0 

J/'0, r? mass etc requirements 

2.6 

Luminosity measurement 

1.4 

MC generator 

1.0 

Trigger simulation 

2.0 

Intermediate decay branching fractions 

1.6 

Signal MC simulation statistics 

0.2 

Fit uncertainty 

64 

Total 

65 


In Refs. [il[il, the partial width of 0(45) —>■ e+e 
was estimated, i.e., F(0(45) —>■ e+e~) = 0.63 keV 15| 
and r(0(45) —>■ e+e“) = 0.66 keV [l 6 [. If we take the 
lower theoretical calculation r(0(45) —>■ e+e“) = 0.63 
keV [T^ . we can obtain the conservative upper limit on 
branching fraction 5(0(45) —>■ 77 J/ 0 ) at 90% G.L., i.e., 
5(0(45) ^77J/0)< 1.3%, 

which does not contradict the theoretical prediction of 

1.9x 10-® il. 
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Table 2. Results of the fits to the e'''e~ cross sections using three coherent resonances: i/)(4040), ^(4160) 

and V'(45'). The errors are statistical only. M, T, and are the mass (in MeV/c^), total width (in MeV), 

product of the branching fraction to and the e"*"e“ partial width (in eV), respectively. <j)i is the relative 

phase between the '!(>(4040) and '!/:(4160) (in degrees) and <f )2 is the relative phase between the ?/)(4160) and 
(in degrees). 



Solution I 

Solution II 

Solution 111 

Solution IV 

-M.|/)(4040) 


4039 (fixed) 


r^(4040) 


80 (fixed) 


S(0(4O4O)^r;J/0)r^|t-"°^ 

6.1 ± 1.1 

11.8 ± 1.4 

12.2 ± 1.5 

6.3 ± 1.1 

-^i/>(4160) 


4153(fixed) 


rVj(4160) 


103(fixed) 


S(0(416O)^7;J/0)rff“°^ 

6.3 ± 1.6 

16.7 ± 2.2 

20.4 ± 2.4 

7.7 ± 1.7 

M^^4S) 


4230(fixed) 


^ip(4,S) 


38 (fixed) 


S(0(45)^t7J/0)rfit- 

0.8 ± 0.6 

0.9 ± 0.7 

4.5 ± 1.3 

4.0 ± 1.2 

01 

320 ± 12 

258 ± 6 

262 ± 5 

324 ± 12 

02 

171 ± 16 

117 ± 17 

142 ± 8 

197 ± 12 


4 Summary 

In summary, we search for the '0(45') state in the 
^ r]J/'tljK^ and e+e“ —>■ ryJ/'f/' processes based 
on the Belle measurements with the assumed mass 
M = (4230 ± 8 ) MeV/c^ and width T = (38 ± 12) 
MeV. The riJ/ij) mass spectrum from B decays and 
the cross sections of e+e“ —>■ r]J/tp are fitted with the 
V’(45) resonance included for the first time. No signif¬ 
icant signal is observed, and 90% C.L. upper limits of 
B{B^ 0(45))-B(0(45) ^ r]J/il;K^) < 6.8 x 10"® and 


B{ip{4S) — 7 > < 8.2 eV are obtained. With 

the r('0(45) —>■ e+e“) = 0.63 keV [l^ as input, we have 
B{'tp{4:S) —>■ < 1.3%, which is almost one order 

higher in magnitude than the theoretical prediction Q. 

The expected integrated luminosity at the Bellell ex¬ 
periment is 50 ab“^ in 2024, which is about 50 times 
the current total integrated luminosity at Belle. With 
this huge data sample, the expected upper limit on 
^(■0(45) —>■ rjJ/'ip) will be 1.9 x 10“® if it scales as 1/VL, 
where L is the integrated luminosity, and can therefore 
reach the theoretical prediction level. 
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